Scaling of elliptic flow in heavy ion collisions 
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Abstract. 

The common interpretation of V2 in heavy ion collisions is that it is produced by hydrodynamic 
flow at low transverse momentum and by parton energy loss at high transverse momentum. In this 
talk we discuss this interpretation in view of the dependence of v% with energy, rapidity and system 
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I size, and show that it might not be trivial to reconcile these models with the relatively simple scaling 

found in experiment 
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INTRODUCTION 



Two experimental results of heavy-ion collisions have been subject to many theoretical 
and phenomenological investigations [1]: One is the observation of a significant sup- 
pression of high- pr particles, "jet quenching", the other one is the observation of an 
azimuthal dependence of the particle spectra on the reaction plane §o n at both high and 
low momenta, the "elliptic flow". The elliptic flow, vj, is parametrized as the second 
Fourier component of the transverse momentum distribution of the produced particles 
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The interpretation of the first finding is generally thought to be that the matter produced 
in heavy-ion collisions is "opaque", with a large energy loss per unit length of fast 
particles [2, 3]. The second finding has been interpreted in terms of the "perfect fluid", 
the hypothesis that matter in heavy-ion collisions has an extremely low viscosity [4, 5]. 
Hence, initial anisotropics in configuration- space density of the collision area will be 
transformed into anisotropics in the collective flow of matter. 

It is important to emphasize that V2 is present at all values of momentum, but has 
different origins at different momenta if the consensus outlined above is correct. Elliptic 
flow at low-/>r is due to the hydrodynamic evolution of the system while elliptic flow at 
high-pr is thought to be due to opacity, since partons emitted in the reaction plane loose 
less energy than partons emitted perpendicular to the reaction plane due to the shorter 
distance traveled. Both are thought to be a dynamical response of the primary asymmetry 
present in heavy-ion collisions. While the scale delimiting these two regimes is assumed 
to be the average pj of the system, (pr) ~ ^(0.5 — 1) GeV (with the tomographic 
regime actually appearing at — 4) x (pr)), the way these two mechanisms combine 



at realistic pj is not entirely clear. In addition, it might be possible that some o/both 
elliptic flow and jet suppression are not generated due to medium but due to initial 
state effects. Strong color dynamics at the parton saturation scale (the "Color Glass 
Condensate") has been shown to exhibit some jet suppression and elliptic flow [6, 7, 8]. 

Phenomenologically distinguishing between different models, even at a given energy, 
is not so trivial because every model has quite a few undetermined fit parameters. Hence, 
for instance, it is not as yet clear whether jet energy loss proceeds by weakly coupled 
or strongly coupled [9, 10] jet-medium dynamics, and we are far from understanding at 
what energy, if ever, do these effects significantly change. 

One important experimental finding which can be used to clarify these questions is 
the discovery of a scaling in elliptic flow across different energies, system sizes and 
centralities, when the data is plotted against transverse momentum p T rapidity ^pseudo- 
rapidity r\ and transverse multiplicity density (l/S)(dN/dy), and eccentricity £. Some 
experimental observations which help us define this scaling are [11]: 

• The dependence of V2/e at mid-rapidity on only the transverse multiplicity density 
(l/S)(dN/dy) across all available energies, system sizes and centralities [12] 

• The "limiting fragmentation" of V2 in rapidity [13] 

• The approximate independence of V2(pr), in a given centrality class, on energy 
[14], system size [15] and rapidity [16] 

Note that the transverse area S and the multiplicity e are theoretical parameters, necessi- 
tating either a Glauber or a Color-Glass event-by-event Montecarlo simulation. Further 
in this work we shall see how to cast some of this scaling in a purely experimental form. 
One way to parametrize all this experimental data, at all energies and rapidities, is 

V2(pt) = £(b,A)F(p T ) , (v2) = e(b,A) J dprF(p T )g (pt, (PT) y ^,b,^ ) (2) 

Here, £ is the eccentricity (dependent, in a Glauber parametrization, somewhat weakly 
on energy, strongly on system size and centrality), g(pr) is the distribution in transverse 
momentum, approximately characterized by one parameter (the average momentum 
(pr) or equivalently, the slope T), which in turn seem to depend, across rapidity y, center 
of mass energy ^ and centrality on just the initial density, in the Bjorken formula . 
F(pr) seems to be a universal function, independent of both energy and eccentricity. 

These are purely experimental statements, with no theoretical overlay, restating the 
results [12, 13, 14, 15, 16] in mathematical form. As such, they are "as good as the 
error bar", and a thorough scan in energy, system size and rapidity might in the com- 
ing years discover violations (some violation of the pj dependence can be seen at low 
Pt in [14]). Taking all this as an established fact, however, is an extraordinarily strong 
constraint, since "typically", for a complicated dynamical model (as non-linear hydro- 
dynamics overlayed with jets inevitably is), generally does not factorize in any way, V2 
is simply V2(xi) (where X{ = {y/s,b,A,y,p T }) and any element of <9 2 V2/<?lnx,(9lnxj is 
non-negligible 

The first scaling, but not the universality of F(pj), was predicted in [17] under certain 
assumptions (no phase transition) within a weakly coupled model (the Knudsen number 
~ 1, so one interaction per degree of freedom per lifetime). Either of the scalings have 



not been thoroughly explored in either hydrodynamics or tomography. In this talk, we 
shall give some qualitative limits to the applicability of each scaling. 



SOME COMMENTS ON SCALING IN HYDRODYNAMICS 

It has long been pointed out, both by heuristic arguments [18, 19] and explicit simula- 
tions [20] that the patterns above pose a problem for the hydrodynamic interpretation of 
V2. Close to the hydrodynamic limit, one expects that V2 is 

• Approximately e since V2(e = 0) = and e small and dimensionless 

• Approximately oc c s (T) since V2(c s = 0) = 0, and the dimensionless c s tracks the 
equation of state 

• V2 is maximum for ideal hydro. Since the Knudsen number Kn, quantifying the 
ratio of the mean free path to the system size, is small and dimensionless, V2 ~ 
yideal ^ 1 _ j^ n ^ j n turn ^ ^ Knudsen number is related to the viscosity over entropy 
density r\/s as well as the system size R,Kn~r\/ (sTR) 

• Finally, v l % eal is a highly non-linear function of the lifetime %uf e , v l 2 eal ~ 
V2(T/ ( y e /To — > oo) x /(t^/to) , which can be numerically shown to be mono- 
tonically saturating, ~ /((pr))tanh(...) in a Cooper-Frye [22, 23] freezeout. T/,y e 
is in turn related to the freezeout temperature and energy density Tf,ef. 

• For p/K <C 1 and isothermal freeze-out, Tuf e / To ~ (eo / e /) 4a > with j | h - Jorken < oc < 

Mhubble depending on how "three dimensional" is the flow. This relation becomes 
more complicated, but qualitatively similar, for systems at high chemical potential. 

In summary, elliptic flow in the hydrodynamic limit should scale as 

*~cj(*™)(l-0(l)±±-) (3) 
e J \Tfx R 2 dyj\ V ' sTR] 

It is clear that only & (Kn) terms mix intensive quantities such as the energy density 
e with extensive ones such as the size R. G (Kn°) "ideal" terms, except for the initial 
time To, depend purely on intensive quantities, giving rise to scaling between systems 
of different sizes. As we will see, this is not true for high pj. For low pj, no scaling 
violation is seen in experimental data, giving a bound for Knudsen number compatible 
with 7]/s = [21], albeit with large error bars (which still need to be computed), at all 
energies. Moreover, the lack of scaling of To is troubling since, by causality, it can locally 
only depend on energy and the local intensive parameters. Just by dimensional analysis, 
it is difficult to see how It can be constant w.r.t. energy, since To\A ~ <^(1)- Landau 
hydrodynamics would imply To ~ 1 / while a CGC type initial condition would most 
likely give a logarithmic dependence since To ~ Qj l . Either, however, would lead to 
unobserved systematic scaling violations. The only way to avoid these is to assume To 
is of the order of the local mean free path at equilibrium, and hence gets reabsorbed as a 
function of the entropy density, (~ (l/S)(dN/dy))~ 1 ^ 3 for an ideal Equation of state. 



Additionally, the Cooper-Frye formula [22, 23] leads to a non-universal F(pj). To 
show this, it is sufficient to expand it azimuthally in eccentricity [22, 23] 
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As long as 8vt,A§p ~ es°, V2(pr) is independent of -y/s. This will be true in the limit 
where the hydrodynamic phase is "long", zy/zb ^> 1, but will not be the case [20, 5] 
if the duration of the hydrodynamic phase < eR, as is the case at lower energies. The 
introduction of an "iso-knudsen freezeout" rather than an isothermal one, a physically 
reasonable scenario explored in [27], should further break this scaling. The reasons for 
this behavior go all the way to the qualitative description of how v 2 behaves in hy- 
drodynamics: V2(pr) and pj integrated (v 2 ) scale differently, because in hydro Fourier 
components of the transverse flow vj(r, 0) depend on lifetime Xuf e differently: 

V2(pr) depends only on the 2nd Fourier component of vj(r, 0) ~ tanh(i/,y e / To) 
(v 2 ) depends on both the 0th ((pr) ~ { T Ufe/ lb) ) and 2nd component. 

Given these, making v 2 (/?r) independent of energy but varying (pr) strongly at all 
energies in unnatural in hydrodynamic models. Detailed simulations including chemical 
potential, however, are needed to determine when does Xuf e become "short", and more 
experimental data might yield a breakdown of V2(pr) scaling at lower energies. 

A further consideration is in order regarding the breakdown in scaling in particle 
species [14]. v 1 2 {pt) = x,K,p,A, ...) does not scale the same way by particle species 
as it does for all particles: different particle species v' 2 (pt) are different at different 
energies, but the differences cancel out when total v 2 ~ Li v 2 (pr)"i(pr) / (Li^i(PT)) 
is considered(n ; is the particle species abundance). In a Cooper-Frye freezeout [22] 
there is no reason for such a cancellation between flow and hadrochemistry to happen. 
Coalescence models, while they will also break coalescence scaling at lower energies 
[24, 26], also do not predict such behavior. 

We note, as a speculative suggestion, the fact that structure functions /(x, Q 2 ) and 
fragmentation functions D q ^i(z,Q 2 ) naturally follow the scaling suggested by both 
the overlap of v 2 (pr) and its breakdown by particle species, since both of these de- 
pend weakly on momentum exchange Q 2 but strongly on the rescaled variables x = 
pP arton /E,z = p hadron I pP arton [25]. In structure functions, x is absorbed into the longi- 
tudinal component of momentum, with Q 2 ~ p\. In fragmentation functions, z ~ pj but 
unitarity protects the effect of fragmentations on all hadrons, / zdzD q ^i(z,Q 2 ) = 1. 
Together, these lead to a particle species dependent v 1 2 (pt) with y/s, but a much weaker 
dependence of total v 2 (pr)- The suggestion that perhaps v 2 (pr) is not a reflection of flow 
at all, but of non-perturbative QCD response of geometry needs much more theoretical 
development, and has to contend with the difficulty of pQCD appearing at y/s = 7.7 
GeV (a comparatively low energy, although much above the "minijet" 1 GeV scale). If 
the overlap of high pj with low pj v 2 , examined in the next section, persists, scenarios 
like this might however need further development. 



HIGHER p T v 2 : SCALING IN THE TOMOGRAPHIC REGIME 



While the scaling we discovered can be, to a certain extent, understood in the hydrody- 
namic regime, the tomographic regime is widely expected to break it. This is because 
in hydrodynamics flow, and hence /^-correlations, is generated by density gradients. 
In tomography ,it is generated by path length variations. Hence, the role of "size", {R), 
which typically depends on system size as ~ A 1 / 3 and is weakly dependent on energy, is 
very different in the tomographic regime w.r.t. the hydro regime. As we saw in the previ- 
ous section, "extensive" factors (RT) (sizex temperature) in the hydrodynamic regime 
are suppressed by 6 (Kn), and hence vanish in the ideal hydrodynamic limit. In the to- 
mographic regime, v 2 ~ ef((RT)) and the {RT) dependence is not suppressed by any 
small parameter. The exact form of the function depends on the details of the model, 
but, as we shall see, it generally does not drop out when {R) and {T) are varied sepa- 
rately by scanning in centrality, system size and y/s. This is true in "standard" jet energy 
loss calculations, both weakly [2, 3] and strongly coupled [9, 10], where opacity is a 
smooth function of the entropy density. It is even more true in models, such as [28], 
where it strongly depends on the distance from the deconfinement transition. Changes in 
the opacity parameter (parametrized by K in Eq. 6) suggested in [32], as well as changes 
in the structure function moments in y/s, are also expected to make the scaling worse. 

Of course, the transition between "hydro" and tomography is a smooth superseding 
rather than a "turning on/off". Indeed, hydrodynamics and tomography can be defined in 
terms of the Knudsen number in momentum space: assuming the scattering cross section 
a depends on the exchanged momentum Q as o ~ l/<2 2 , and assuming momentum 

is much higher than temperature, the Knudsen number becomes ~ ^ Hence, 
the tomographic regime starts dominating when Kn ~ 1. For different energies and 
transverse densities the critical pr can be easily shown to be 

p'rTiV^MAi) (S 2 dN l /dy \ a 
Pt2° [V^MM) ~ \Si dN 2 /dy) 

where CO = 1/2 for collisional-dominated equilibration but could increase to 3/2 [29] if 
radiative processes become important in the hydrodynamic regime. 

The advantage of characterizing as "hard" hadrons with pj > pj mo is that this defini- 
tion is independent of details such as the dynamics of production and fragmentation of 
fast hadrons (we do not have to call them "jets", which is problematic at low energies). 
Experimentally, the fact that p^ mo decreases with decreasing N part , y/s is advantageous. 

v 2 in the tomographic regime of course continues to be a good observable. In particu- 
lar, it is independent of effects like Cronin effect and jet reconstruction, which become 
problematic at low energies: At lower energies Raa diverges, due to kinematics, at any 
medium opacity. Since kinematic effects, by themselves, do not depend between the 
hadron-reaction plane angle, however, v 2 does need opacity to be generated. This makes 
comparing v 2 (pr > Pj mo ) at different energies and system sizes, and looking for scaling 
violations, an optimal probe of changes in opacity with temperature. In the next section 
we use the ABC model [31, 32] to put some limits on this scaling in a class of models 
where v 2 (pt ^ (pt)) is generated by tomography alone. We shall make the case that an 
experimental investigation of this scaling is a useful test for energy loss models. 



A theoretical investigation using the ABC Model 



Since the reaction plane was determined at low pr, one can now eliminate the theoret- 
ical input e, S from our observables by concentrating on V2(pr) / {vi}- The denominator 
is the momentum-integrated V2, which should include all eccentricity information. If V2 
depends purely on gradients, of course, this ratio should be independent of both energies 
and system sizes. The size dependence, however, should lead to a break in the scaling. 

The ABC model is a simple parametrization which describes the energy loss of a 
"fast" particle (pr/T <C 1) in traveling "large" medium (1 / (Ti) <C 1, T is the propaga- 
tion time). If the parton is light and on-shell, energy loss models should give 

ax \pt I t 

a, b, c are a nice phenomenological way of keeping track of every jet energy loss model 
in certain limits. In a collisional dominated parton cascade a = 1 , c = [30], in a radiative 
dilute plasma ("Bethe-Heitler regime") a = 1/3, c = 0, in a dense plasma (LPM regime) 
a = l,c = 1) while in a "falling string" AdS/CFT scenario [9, 10] a = 1 /3,c > 2. 

We can now expand in "empirical" parameters (^) , £, and get 

v 2 {pf h ) ~ B«L°*W) P F {aAc) adn ' dy . (7) 

with non-trivial exponents a ; , calculated in each model in [33]. The size parameter 
L ~ min ((/?) , zy — To) , sampling the lifetime for some systems and the transverse size 
(R) for others. We perform a numerical investigation with a background given by a 
longitudinally expanding ellipsoid, fitted to global multiplicity and system size 

1 /3 

T(rr, 0, T, To) = 7b (^) (r -7? (1 + ecos (20))) (8) 

Where R, e are scanned scanned across radii of Cu,Au,Pb, To are chosen according to 
the assumptions described in section To = \fm,^fs 1 ,T Q j -1 , and Tq is adjusted to re- 
produce multiplicity and all energies. We then obtain V2(pr)/ (V2), the latter given by 
the experimental parametrization [vi) /e = 4(10~ 3 ) (1 /S)(dN/dy). The results, shown 
in Fig. 1 and 2 show that scaling fails, by a similar magnitude, in all models. Currently, 
experimental data does not allow to make definite conclusions but, as shown in Fig. 3 
[12, 34, 15], modulo rather big error bars, the only scaling violation is seen at interme- 
diate regions. Scaling holds across different centralities up to well above pr ^ 20 GeV, 
and seems to break up only at p T ~ 40 GeV at the LHC. While a systematic shift of the 
center is seen comparing RHIC and LHC, the error bars are way too big to attach any 
meaning to this conclusion. This shift is not seen up to p T ~ 3.5, where, as noted in 
[15], this scaling holds across RHIC Cu+Cu and Au+Au. 

If such a result survives when the error bars decrease, and persists across Au+Au vs 
Cu-Cu at RHIC and Pb-Pb vs Ar-Ar at the LHC, it would be a very profound statement 
given Figs 1 and 2. It would make it inevitable that non-tomo graphic effects play leading 
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FIGURE 1. Scaling plots assuming z = 1, indicative of a pQCD dense plasma (LPM limit) energy loss, 
for RHIC Cu+Cu,RHIC Au+Au and LHC Pb+Pb 
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FIGURE 2. Same as the previous figure, but z = 2, indicative of AdS/CFT energy loss 



roles in high pj vi- While such effects have been suggested , for both the initial stage [8] 
and final fragmentation [35], it is generally not expected that they play a large role at high 
Pt- A continuation of scaling might force us to revise such expectations, since causality 
makes it inevitable that short-lived processes scale as gradients (such as e) rather than 

extensive quantities such as (R a T^ \ We therefore eagerly expect further scaling studies 

at high pr V2 to test our expectation that vi{pt 3> (pr)) is generated tomographically, 
and to use the scaling violation as a constraint on jet energy loss models. 
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